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Abstract 

The space-time conservation element and solu- 
tion element (CE/SE) method is used to solve 
the conservation law form of the compressible 
axisymmetric Navier-Stokes equations. The 
equations are time marched to predict the un- 
steady flow and the near-field screech tone 
noise issuing from an underexpanded circular 
jet. The CE/SE method uses an unstructured 
grid based data structure. The unstructured 
grids for these calculations are generated based 
on the method of Delaunay triangulation. 

The purpose of this paper is to show that an 
acoustics solution with a feedback loop can be 
obtained using truly unstructured grid technol- 
ogy. Numerical results are presented for two 
different nozzle geometries. The first is con- 
sidered to have a thin nozzle lip and the second 
has a thick nozzle lip. 

Comparisons with available experimental data 
are shown for flows corresponding to several 
different jet Mach numbers. Generally good 
agreement is obtained in terms of flow physics, 
screech tone frequency, and sound pressure 
level. 

1 Introduction 

It is desirable to develop a general computer code that 
can predict the complex flow structure around compli- 
cated geometries and the associated aeroacoustics. Many 
computational fluid dynamics methods have been devel- 
oped that take advantage of an inherent grid structure 
which permits flow solutions to be obtained efficiently. 
For computational aeroacoustics (CAA), another obvi- 
ous advantage of the structured grid is its well-controlled 
grid quality (size and shape) that helps to maintain the 


high quality of the computed acoustic waves. This struc- 
ture, which makes the solver so efficient, often makes 
it difficult to generate grids about complex geometries 
without an additional complication of blocked or overset 
grid techniques. 

An unstructured grid flow solver can alleviate many of 
the problems associated with structured grids. Unlike a 
structured grid, the unstructured grid neighborhood must 
be provided to the solver explicitly. The simplest geo- 
metric figure that can cover a two-dimensional domain is 
the triangle. The advantage of the triangle is that it can 
easily be used to generate a grid over an arbitrary domain 
without overlapping. It also has the advantage of cluster- 
ing in high flow gradient regions without the concern of 
the surrounding cell structure. The main disadvantage of 
the unstructured grid is the additional memory require- 
ments in the flow solver to store the grid. More of the 
burden is placed on the computer code, thus the user in- 
put is reduced both during grid generation and flow so- 
lution such that more configurations can be tested in a 
shorter period of time. For CAA, the grid quality of an 
unstructured grid is less controllable and it may deterio- 
rate the acoustic wave computations. The situation may 
be even more critical if a complicated acoustic feedback 
loop is present, which is the problem type considered in 
this paper. 

The problems under consideration in this paper are cir- 
cular underexpanded supersonic jets [1-8]. These jets ra- 
diate a combination of mixing noise, broadband shock- 
associated noise, and under certain conditions screech 
tones. Mixing noise from jets can be directly associ- 
ated with large-scale structures, or instability waves in 
the shear layer. The broadband shock-associated noise 
and screech tones are associated with the interaction of 
these instability waves with shock cell structure in the 
jet core. The screech tones arise due to a feedback loop, 
i.e., part of the acoustic waves generated by the instabil- 
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ity wave/shock-cell interaction propagate upstream and 
re-generate the instability waves at, or in the vicinity of, 
the nozzle lip. The desire to control screech noise is 
not only due to general noise reduction concerns. It can 
also cause damage to the structure due to sonic fatigue. 
The feedback loop is sensitive to small changes such as 
lip thickness and other phenomena. Most of this under- 
standing comes from experimental observations. Numer- 
ical methods that can predict screech tones will help in 
the understanding and the eventual control of this noise 
source. 

A direct numerical simulation was computed by Shen 
and Tam [9,10] for screech tones issuing from a circular 
jet, based on the well-known DRP scheme, with a struc- 
tured grid. A large computational domain of 35D long 
in the streamwise direction and 17 D in radial direction 
was used in their computation, with D being the nozzle 
diameter. 

More recently the space-time conservation element 
and solution element (CE/SE) method with a triangulated 
structured grid [11-14] has been used to predict the near- 
field screech tones from underexpanded supersonic jets 
[15]. The near-field results compared well with the ex- 
perimental data, [4-8]. This work also demonstrated the 
capability of the CE/SE nonreflecting boundary condi- 
tion. This allows the computational domain to be more 
compact reducing the overall mesh memory size require- 
ments. The present work is a companion work of [15]. 
Here a truly unstructured grid, generated by a point in- 
sertion Delaunay method rather than a triangulated struc- 
tured grid, is employed. This is done with the idea that 
more complex geometries will require a more flexible 
grid generation technology. The purpose of this paper is 
to show that an acoustics solution with a feedback loop 
can be obtained using unstructured grid technology. 

The governing equations are given in Section 2. Here 
the nondimensionalized unsteady axisymmetric Navier- 
Stokes equations in conservation law form are presented. 
Unstructured grid generation is discussed in Section 3. 
The requirements of the grid for the flow solver will be 
explained. Also the flow code was executed on a parallel 
computer architecture, so the domain decomposition and 
message passing libraries are briefly reviewed. The un- 
structured based axisymmetric CE/SE method used here 
is briefly discussed in Section 4. Section 5 will provide 
the problem description with the initial and boundary 
condition requirements for the particular screech prob- 
lems examined in this paper. In Section 6 numerical re- 
sults are presented and and compared to previous numer- 
ical findings using structured grids and available exper- 
imental findings [4-8]. Conclusioning remarks will be 
made in Section 7. 


2 Governing Equations 

The unsteady axisymmetric Navier- Stokes equations 
are used to model the viscous problems in this study, 
since experimentally for the thin lip nozzle [4-8] with 
Mj < 1.19 the flow is axisymmetric. In nondimension- 
alized differential form the conservation laws are written 
as: 


Ut + F x + G y -i Q, (1) 

where x , y, and t are the streamwise and radial coordi- 
nates and time, respectively. The conservative flow vari- 
able vector, U, and the flux vectors in the streamwise 
direction, F , and radial direction, G, are given by: 


U 


with 



Us 

\U 4 



G = 



Ui = p, U 2 = pu , U 3 = pv, 
04 =p/( 7 - 1) + p(u 2 + v 2 )/ 2 . 


Here p,u,v, p , and 7 are the density, streamwise velocity 
component, radial velocity component, static pressure, 
and constant specific heat ratio, respectively. The flux 
vectors are further split into inviscid and viscous fluxes: 


F = F i -F v ,GmG i -G v , 


where the inviscid fluxes are given as: 


Fi 1 = U 2 , 

F a = (7 - l)tf 4 + [(3 - 7)0| - (7 - 1)0?] /2U U 
F i3 = 0203/01, 

Fi 4 = 702 04/01 - (7 - 1)02 [0 2 2 + 03 2 ] /201, 
Gn = 03 , Ga = 0203 / 01 , 

Gi 3 = (7 - 1)04 + [(3 - 7)0 3 2 - (7 - 1)0 2 2 ] /201, 
G i4 = 70304/01 - (7 - 1)03 [0f + 0 3 2 ] /201 2 , 

and the viscous fluxes are: 


F v 1 = 0 , F v2 = p(2u x - §V-V), 

Fy3 = P j (v x + r U j y)i 

F v 4 = /i [2 uu x # (u y + v x )v - §(V*V)tH- 
7 d U 4 u 2 + v 2 

P^dy^Th 2 ^ 

G v 1 = 0 , G v 2 = p(v x + u y ), 

G v 3 = p(2v y — |V*V), 
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G v 4 = /i[2 VVy + (% Hb Vaj)^ - §(V-V)u+ 

7 9 U± u 2 +v 2 
P^dy^l h 2 ^ 

where u,v,u x ,u y ,v x ,v y are respectively the x and ?/ 
flow velocity components and their derivatives, which 
can be written in terms of the conservative variables 
Ui,U 2 ,U 3 and U±. The Prandtl number and viscosity 
are Pr and /i, respectively. The velocity divergence is 


V-V = u x + v y + v/y. 


The source term, Q, appearing on the right hand side 
of the equation, Eq. (1), is a result of the axisymmetric 
formulation and is given as: 


Q 



where 


Q i = —U^/y, Q2 = -^U 2 Us/Uiy, 

Q3 = -Ui/lhy, Q 4 = -G 4 /y. 

Considering (x,y,t) as coordinates of a three- 
dimensional Euclidean space, £ 3 , the Gauss divergence 
theorem can be used to give the following equivalent in- 
tegral form of the conservation laws, Eq. (1): 




m = l,2,3,4, (2) 


where S(V) denotes the surface of the volume V in £3 
and H m = (F m , Gm , U m ) represents the total flux leav- 
ing the boundary of space-time region. The integral form 
of the conservation laws allows the capture of shocks 
and other flow discontinuities naturally. This integral 
form of the equations is solved numerically by the space- 
time conservation element and solution element method 
which is described in Section 4. 


3 Grid Generation 

An unstructured grid is used to show that aeroacous- 
tics with a feedback loop can be predicted in complex 
domain shapes. A triangular shaped cell will be used 
as the geometric shape to cover a two-dimensional com- 
putational domain. The use of an unstructured grid re- 
moves the available implicit reference of the local com- 
putational neighborhood by incrementing an index as in 
a structured grid. Instead, the neighborhood is provided 
explicitly through a connectivity matrix. This connectiv- 
ity matrix usually contains cell based information as well 
as edge based information. 



Figure 1: Triangulation of computational domain with 
enlargement showing nozzle lip region. 


The use of a triangular unstructured grid formulation 
has some distinct advantages over a structured grid. One 
advantage to using triangular cells is that with them it is 
easy to generate grids about complex geometries. This 
reduces the amount of time required to generate a suit- 
able grid. Also grid adaptation can be done locally with- 
out adding unnecessary cells to other regions of the do- 
main. 

The method used in this work follows a path similar to 
that of Holmes and Snyder [16] to triangulate a region. 
An example of an unstructured grid for the thin lip nozzle 
geometry generated in the above manner is shown in Fig. 
1. The computation is axisymmetric, so only the upper 
half of the domain needs to be considered. An enlarge- 
ment of the nozzle lip region gives an indication of the 
grid resolution used in these aeroacoustic feedback loop 
problems. 

The flow code requirements dictate the type of output 
that the grid generation scheme must provide. A connec- 
tivity array must be generated for an unstructured grid 
so that a cell neighborhood is completely defined for the 
flow code. The code in the current work is based on a 
cell centered scheme. Here the triangle itself is the con- 
trol volume used in the finite volume formulation. 

Connectivity is determined by cell nodes, cell edges, 
and edge cells. Cell nodes are the nodes at the vertices 
of a triangle. The cell edges are the edges of a triangle 
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that connect those nodes. The edge cells are the cells 
that lie on either side of a given edge. Typically the x , 
y coordinates of the cell nodes are the only grid floating 
point numbers required as input by a flow code. The grid 
connectivity is defined through integer arrays. 

Boundary condition information must be defined ex- 
plicitly. Ghost cells are used across boundaries such that 
specific boundary conditions can be applied. A boundary 
condition such as solid wall, inlet, exit, or any other user 
defined condition can be prescribed along any cell edge 
through the ghost cell. 

For aeroacoustic problems, large numbers of cells are 
needed to resolve the unsteady waves. The grids are 
made more manageable through the use of parallel com- 
puter architecture. The publically available code METIS 
[17] is used to decompose the unstructured mesh into a 
predetermined number of subdomains so that it can take 
advantage of computers with multiple processors. The 
mesh partitioning code then load balances, minimizes the 
communication between processors, and minimizes the 
maximum number of subdomains adjacent to any other 
subdomain. Each processor then operates on a subset 
of the larger domain. It is necessary to communicate the 
subdomain boundary information to the neighboring pro- 
cessor so that the numerical algorithm can perform com- 
putations that require data from its immediate neighbor. 
The message passing libraries provided by the standard- 
ized package MPI [18] is used for point-to-point com- 
munication between processors. For unstructured grid 
based solvers a list must be maintained that provides the 
boundary cells so that the appropriate data may be sent 
and received by a processor. 

A windowing technique is used to pack points in re- 
gions of the computational domain where more resolu- 
tion is desired. To a limited extent the smoothness of the 
grid is sacrificed for an increase in resolution using this 
technique. A grid generated by this method will be re- 
ferred to as windowed Delaunay point insertion (WDPI). 

The unstructured grid generator Gridgen [19] is also 
used to generate a grid in the same domain for compari- 
son purposes. The stretching techniques in this grid gen- 
erator results in a much smoother grid than that provided 
by the windowing technique. 

4 Discretization Technique 

The space-time conservation element and solution el- 
ement method is a finite volume method where the el- 
ements include space and time coordinate dimensions. 
The original two-dimensional CE/SE method was de- 
rived for triangular space-time grids. To date the grids 
used in computing the aeroacoustics from jet flows us- 
ing this method have been triangulated structured grids. 
Here the method is applied to more general unstructured 
grids generated by the Delaunay method described in the 


t 


Figure 2: Conservation elements used for updating solu- 
tion. 

previous section. In this section a general review of the 
CE/SE method is given as it is applied to jet noise calcu- 
lations. 

There are three conservation elements (CE’s)that are 
associated with a given triangular space- time element, 
Fig. 2. The three conservation elements are shared ge- 
ometrically by the three adjacent triangular space-time 
elements. The geometric center of the triangular faces 
adjacent to face A ABC are shown as points D , E , and 
F. Each conservation element has six faces that form 
a quadrilateral cylinder in space-time. The bottom and 
top faces are formed by connecting the centroids of the 
triangular bottom and top faces (of the space-time trian- 
gular element) to the nodes that make up that face, re- 
spectively. For a given space-time triangular element, 
this construction produces three quadrilateral faces on 
each of two spatial planes (i.e., A , D, B , O ). These are 
nonoverlapping elements. For a given space-time trian- 
gular cylinder, the union of the three conservation el- 
ements forms a hexagonal cylinder of volume V (the 
lower face of this hexagonal cylinder is composed of 
nodal points A,D,B,E,C,F). The conservation law 
Eq. (2) is applied to this hexagonal cylinder. The result- 
ing solution is associated with the center of the hexagon 
at the new time level (not necessarily coincident with 

O'). 

The solution elements (SE) used here are the inter- 
faces of the conservation elements that form the hexago- 
nal cylinder, Fig. 2. 

The integral equation, Eq. (2), will provide a solution 
at the geometric center of the top hexagonal face. This 
does not necessarily coincide with the triangle cell center 
(the solution element) location, E. Discretization of the 
integral conservation laws, Eq. (2) leads to 

Y J H* m -^ = V{Q m ) n +\ (3) 

s(Y) 
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NRBC 


for m = 1, 2, 3, 4, where the space-time flux at the sur- 
face centers of S(V), H ^ = (T^, G^, U^), can be ob- 
tained. The source term on the right-hand side of Eq. (3) 
is computed at the center of the hexagon A! D' B' E'C' F' 
at the new time level. 

To evaluate the resulting system of algebraic equa- 
tions, Eq. (3), the values must be computed at the solu- 
tion element faces (actually the face centers). These face 
centered values are approximated by the linear Taylor se- 
ries expansions. The unknowns in the integral equation, 
Eq. (3), are located at the cell center of the hexagon as- 
sociated with O'. Expanding the integral equation for a 
given space-time element results in 12 scalar equations 
for the 12 unknowns at O ' . Each of the three conser- 
vation elements contributes four scalar equations. More 
details of the solution method can be found in [14]. 

The simple non-reflecting boundary conditions de- 
scribed previously [20,21, 15] still work well with a De- 
launay based unstructured grid. The weighted a — e 
CE/SE scheme is used in this work. To suppress wig- 
gles around shocks the flux derivatives use a weighted 
average similar to that of van Albada [22]. Details are 
provided in [14]. 

The source term evaluation was presented in a previ- 
ous paper [15]. Briefly, the source term Q is a function 
of the unknown U\ therefore, a local Newton iterative 
procedure is needed. The procedure takes about 2-3 iter- 
ations to converge. 

A simplified implicit LES procedure similar to those 
used in [23] is adopted here to account for the strong 
momentum exchange in the shear layer. In the LES, a 
simple Smagorinsky subgrid scale model is used for the 
eddy viscosity: 

A h = {C s A) 2 (2SijSij) 1 ^ 2 , 


where 


1 (SUj 
2[ dx j 


du^ 

dxi 


with A = (AxAy) 1 / 2 and C s = 0.1. Then, /i + fi t 
replaces /i in the actual computation. In practice, the 
LES model was only applied for radii less than 1.5 D and 
4 D downstream of the nozzle exit plane. 


5 Problem Description 

The physical domain of interest for the screech prob- 
lem is shown in Fig. 3. A slice of the axisymmetric do- 
main shows the thin and thick lip nozzle geometries and 
the near-field boundaries where boundary conditions are 
specified in the flow code. The circular jet flow enters 
through the nozzle from the left and expands into the sta- 
tionary ambient atmosphere. The flow condition speci- 
fied at the nozzle exit at the start of the computations is 
choked. 



Figure 3: Geometry of the computational domain. 


Here we are interested in predicting the source of the 
noise rather than its propagation into the far-field. The di- 
ameter, D , of the jet nozzle is chosen as the length scale. 
The density, po> speed of sound, ao, and temperature, To, 
in the ambient flow are taken as scales for the dependent 
variables. 

The Delaunay grid generator is used to cluster the tri- 
angular cells near the nozzle lip to provide good resolu- 
tion for the shock cells and the acoustic feedback loop 
path. Cell stretching is used for these computations. We 
will compare flow solutions computed on a Delaunay 
based unstructured grid, described in Section 3 of this 
paper, with those previously computed on a triangulated 
structured grid [15]. The triangulated structured mesh is 
constructed by dividing each rectangular cell of the struc- 
tured mesh into four triangles. An unstructured connec- 
tivity matrix is then produced. More specifics about the 
numbers of cells and flow code parameters will be pro- 
vided in the the results section of this paper. 

To clearly display the upstream propagating screech 
waves, the computational domain was extended 2D up- 
stream of the nozzle exit. The full computational do- 
main is a circular cylinder of 8D axial length and 6D 
radius. At the nozzle exit, the inflow plane is recessed 
by two cells so as not to numerically restrict or influ- 
ence the feedback loop. The first nozzle geometry def- 
inition follows that given by Ponton et al. [8] where a 
straight nozzle lip of 0.2 D in thickness is adopted. The 
second nozzle geometry has the same computational do- 
main, however the nozzle lip thickness is increased to 
0.625D. The computational domain is modified by ex- 
tending the upstream, downstream, and radial boundaries 
by 2D to determine if the boundary locations has an ef- 
fect on the solutions. However, no discernible effect is 
realized. Only one case was used to test the effect of the 
extended computational domain. 


NASA/TM— 2002-21 1799 


5 



5.1 Initial Conditions 

Initially, the flow in the computation domain is set to am- 
bient flow conditions, i.e., (using nondimensional vari- 
ables) 


while U'j and ( U y )™ are now defined by simple extrap- 
olation from the nearest interior node j', i.e., 

U] = u]r 1/2 (U y )] = (U y )]r 1/2 . 


PO = 1 , 



Uo = 0 , 


vo = 0 . 


5.2 Boundary Conditions 

At the inlet boundary, the conservative flow variables and 
their spatial derivatives are specified to be those of the 
ambient flow, except at the nozzle exit, where an elevated 
pressure is imposed, i.e., the jet is under-expanded, as in 
the physical experiments. By using the ideal gas isen- 
tropic relations, it follows that the nondimensional flow 
variables at the nozzle exit, with M e = 1, are given by 


Pe = 


7(7 + !)Pe 

2T r 


2 + (7 — 1)M? 1 7-1 

7 <ir 1 

= (Jri) ’ Ve = 0, 

where T r is the reservoir (plenum) temperature. We will 
also follow the experimental cold-flow condition where 
the reservoir temperature is set to the ambient tempera- 
ture, i.e., T r = 1. 

At the symmetry axis, i.e., y = 0, a simple reflective 
boundary condition is applied. At the top and outflow 
boundaries, the Type I and Type II CE/SE non-reflecting 
boundary conditions as described in the next subsection 
are imposed, respectively. The no-slip boundary condi- 
tion is applied on the nozzle walls. 



5.3 Non-Reflecting Boundary Conditions 

In the CE/SE scheme, non-reflecting boundary condi- 
tions (NRBC) are constructed to allow fluxes from the 
interior domain of a boundary CE to smoothly exit the 
domain. There are various implementations of the non- 
reflecting boundary condition and in general they have 
proven to be well suited for aeroacoustic problems [15]. 
The following NRBCs are employed in this paper. 

For a grid node (j, n ) lying on the outer radius of the 
domain, the non-reflective boundary condition (Type I) 
requires that 


(U X )] = (Uy)] = 0, 

while U™ is kept fixed at its initial steady boundary 
value. At the downstream boundary, where there are 
substantial gradients in the radial direction, the non- 
reflective boundary condition (Type II) requires that 

(U x )] = 0 , 


These are identical boundary condition specifications 
used in [15] for the results computed on a triangulated 
structured grid. 


6 Results 

The CE/SE method is applied to the underexpanded 
axisymmetric circular jet, with two nozzle lip thick- 
nesses. The thin lip nozzle has a wall thickness of 20% 
of the nozzle jet exit diameter. This case is used to com- 
pare the results obtained on a Delaunay based unstruc- 
tured grid with experimental data [ 8 ] as well as with re- 
sults previously obtained using a triangulated structured 
grid [15]. The flow solver for both grids is identical. 
The flow data is compared at a jet Mach number, Mj, 
of 1.19. Screech frequency and sound pressure levels 
over a range of jet Mach numbers are also compared for 
this thin lip nozzle configuration. The thick lip nozzle 
has a wall thickness of 62.5% of the nozzle jet exit di- 
ameter. This case is compared with expermental results 
[ 8 ] over a range of jet Mach numbers. A large number of 
times steps are required to achieve appropriate accuracy 
in the Fourier analysis of the time series data. Typically 
80K time steps are needed to remove the initial tran- 
sients from the computational domain. Then, another 
330AT time steps are required to record the time series 
data. 

It is important to emphasize that no harmonic forcing 
is imposed in the numerical simulation. The initial im- 
pact of the boundary condition at the nozzle exit stim- 
ulates the jet shear layer and triggers the feedback loop 
that generates the screech waves. 

In the following, the computed jet flow data and acous- 
tic screech data are presented separately, and compar- 
isons are made with available experimental results. The 
thin lip nozzle is examined first and comparisons are 
made with the experimental and the previously obtained 
numerical results. The results of the thick lip nozzle are 
then presented and compared with the experimental re- 
sults. 

6.1 Flow and Shock-Cell Structure for the Thin 
Lip Nozzle 

For the flow data, the CE/SE method is applied to the 
thin lip nozzle case at a Mj of 1.19. The resulting shock 
cell structure for the axisymmetric flow calculation af- 
ter 200AT time steps is shown as a Schlieren (density- 
gradient modulus) contour plot, Fig. 4. The experimen- 
tal [4] shock cell structure is shown in Fig. 5, and the 
numerical results show good agreement with the shock- 
cell structure in the experiment. For example, the shock 
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Thin Lip Nozzle 



Figure 4: Instantaneous snapshot of shock cell structure 
for thin lip nozzle, (numerical Schlieren) Mj = 1.19. 
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Figure 5: Experimental Schlieren plot showing shock Fi S ure 6: Frequency of thin lip nozzle over a range of 
cell structure, Mj = 1.19. 


cell width is about 0.81} in the stream wise direction. In 
both plots, it is observed that the first two shock cells 
appear to be sharp and clear since the shear-layer insta- 
bility wave is too weak at these locations to significantly 
affect the shock cell. However, it interacts strongly with 
the shock cells once it has gained a sufficient amplitude 
through its streamwise growth, which is evident by the 
deformation of the third and fourth shock cells. The de- 
formation and movement of these shock cells are also 
seen in the experiments [4,5]. This is, as pointed out by 
Seiner [2], the mechanism by which the acoustic waves 
that produce the broadband shock-associated noise and 
the screech tones are generated. Furthermore, the fifth 
shock cell and those further downstream are so severely 
deformed that they almost disappear. 

As demonstrated in the previous numerical results 
[15], the screech wave length of 1.6D is predicted in the 
computations. This is in good agreement with the exper- 
imental results (e.g.[ 4]) and the theoretical predictions of 
Seiner [2] and Tam [1]. 

The thin lip nozzle flow solution on the WDPI and 
Gridgen based grids compare well with the solution ob- 
tained previously on the triangulated structured grid [15]. 
The triangulated structured grid used a grid consisting of 
132 K cells. The WDPI grid has 167 K cells. The Grid- 
gen has a 13 6K cell grid. Since the flow code is the 
same for all of the grids presented here, there is no clear 
advantage to using a triangulated structured grid. In fact, 
a truly unstructured grid should be much more effective 
when the CE/SE method is applied to complex geome- 
tries. 


6.2 Screech Frequency and Sound Pressure 
Level 

Because the data is unsteady, many of the same tech- 
niques used to analyze experimental data are also used 
to analyze the numerical results. At a point on the nozzle 
exit lip (x = 2.0, y = 0.642) the time history of pressure 
is stored over a large number of time steps. This data 
is then post processed by Fast Fourier Transform (FFT) 
techniques to obtain spectral information. 

The resulting screech frequencies are plotted over a 
range of jet Mach number in Fig. 6. Two experimental 
data sets are provided by Panda et al. [6] and Ponton 
et al. [8]. The nozzles in the experiments have the same 
diameters and operating Mach number, however they dif- 
fer slightly in shape and thickness. The geometry used in 
the numerical predictions is that of Ponton et al. [8]. The 
CE/SE method is used to compute a solution on the same 
jet geometry but on three different grids. The first series 
of solutions are computed on the triangulated structured 
grid [15]. The second set of numerical data are computed 
on a grid generated by WDPI. This screech frequency 
compares quite well with that in [15]. They slightly un- 
derpredict the experimental data at the lower jet Mach 
numbers but they do predict the same trend. At the higher 
levels of Mj the screech frequencies compare well with 
experiment. 

The third set of solutions are computed on a grid gen- 
erated by Gridgen [19]. The frequencies produced by the 
solution on the Gridgen grid are lower than those pro- 
duced on the the other two grids. At a low jet Mach num- 
ber, Mj , of 1.09 there are no detectable acoustic waves. 
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Figure 7: Sound pressure level of thin lip nozzle over a 
range of Mj . 


The computations at the other jet Mach numbers under- 
predict the A i mode screech frequencies observed in the 
experiment, but they do predict the same trend. The A 2 
mode screech frequencies give a close match to the ex- 
perimental data. 

The sound pressure level (SPL) is shown in Fig. 7. The 
computed results are in general lower than the experi- 
mental data. They are closer to the experimental results 
of Ponton et al [8]. The solutions on the WDPI grid 
compare well with the triangulated structured grid. The 
solutions on the grid generated by Gridgen [19] under- 
predicts the SPL of the A 2 mode significantly. This dif- 
ference could be due to the lack of resolution in that grid 
downstream of the nozzle. 

In general it is felt that it is necessary to cluster cells 
in critical flow areas. The grid generated by WDPI is 
probably better in this case since points are forced into 
the region of the flow where the interaction between the 
shock cell and shear layer is strongest. A careful choice 
of a and time step is necessary to obtain the appropriate 
amount of dissipation. 

6.3 Flow and Shock-Cell Structure for the 
Thick Lip Nozzle 

For the thick lip nozzle case with a Mj = 1.19, the 
computed shock cell structure is shown in Fig. 8. This 
is an instantaneous snapshot of isobars combined with 
a numerical Schlieren. Unlike the thin lip nozzle case, 
the shock cell location which interacts with the instabil- 
ity wave begins between the second and third shock cell. 
This is much closer to the nozzle lip than is observed in 



Figure 8: Instantaneous snapshot of screech wave radia- 
tion for thick lip nozzle, (isobars and Schlieren) Mj = 
1.19. Solution shows a strong subharmonic. 


the thin lip nozzle case. This is probably due to a stronger 
feedback mechanism due to the thick lip. Figure 8 is the 
larger computational domain, alluded to earlier in this 
paper, to test the effect of boundary location. The nozzle 
exit location is not changed, so the pressure data is still 
taken at x = 2.0, y = 0.642. 

The thick lip case shows a strong fluid flow entrain- 
ment which produces a counter rotating vortex just up- 
stream of the nozzle lip, Fig. 8. This entrainment vortex 
is observed for over 500 K time steps. The location of 
the vortex practically does not move. The vortex passes 
cleanly through and breaks the acoustic waves as they 
move upstream. 

A strong subharmonic wave is observed in this calcu- 
lation. It propagates from a region in the flow where the 
pairing of vortices occurs and out through the upstream 
boundary. The subharmonic wave is also split by the en- 
trainment vortex. Due to the two-dimensional axisym- 
metric aspect of the computation, the subharmonic wave 
may be much stronger than the physical one. 

The total number of grid cells for the WDPI grid is 
184 K cells. The grid is clustered in the region of the 
nozzle lip using the windowing technique described ear- 
lier. The Gridgen grid contains 140 K cells. This grid is 
also clustered but with smoother results than that given 
by windowing. Figure 8 has 250 K cells with grid clus- 
tering in the vicinity of the nozzle lip. 


NASA/TM— 2002-21 1799 


For this geometry an a of 0.5 is used for the WDPI 
solution with a time step of 0.0005. Any larger time step 
requires a correspondly larger a to the point where it af- 
fects the acoustic solution. The Gridgen based solution 
still could be run with an a of 1.0 and a time step size of 
0 . 002 . 

6.4 Screech Frequency and Sound Pressure 
Level 

The FFT is used to analyze the unsteady feedback loop 
that is generated between the shock cell structure and 
the nozzle lip. The pressure history is recorded over a 
large number of time steps, Fig. 9. The pressure level 
is recorded every twentieth time step after a time period 
such that the start-up transients have washed out of the 
computational domain. Figure 9 represents the time trace 
of pressure (nondimensional) at the point on the nozzle 
exit lip, x = 2.0, y = 0.642 for a Mj = 1.15 using the 
WDPI grid. This signal is then analyzed by FFT. The 
resulting SPL spectrum is shown in Fig. 10, where only 
the frequency range of 0 to 20000 Hz is displayed. There 
are two SPL spikes in the spectrum that correspond to es- 
timated screech frequencies of 7195 Hz (SPL = 142dB, 
Ai mode) and 9169 Hz (SPL = 136dB, A 2 mode) with 
a binwidth (or bandwidth) of 326 Hz. This indicates that 
there is either two screech tone frequencies occurring at 
the same time or that screech mode switching is present. 
A strong subharmonic is also shown in this case. 

The screech frequencies at the lower jet Mach num- 
bers matched the results of Ponton et al. [8] very well, 
Fig. 11. At higher jet Mach numbers, Mj = 1.15,1.17, 
and 1.19, the screech frequencies are overpredicted. 

The sound pressure levels for the thick lip nozzle in the 
experiment has a higher level for the A 2 mode than for 
the Ai mode, Fig. 12. This is not seen in the numerical 
predictions. 

A major advantage of the grid produced by Gridgen 
is that it requires less input from the user than does the 
WDPI grid. The grid resolution could be augmented con- 
siderably in the feedback region of the flow domain with 
proportionately better results. The smoothness charac- 
teristics of the Gridgen grid may also play an important 
role in getting a better solution. 

The axisymmetric computation may be predicting a 
much larger subharmonic frequency than would be seen 
in an experiment or in a three-dimensional computation. 
This may be due in part to the fact that the flow is being 
artificially confined as axisymmetric. 

7 Concluding Remarks 

In this paper, the unstructured CE/SE Navier- Stokes 
method is applied to underexpanded supersonic axisym- 
metric jets over a range of jet Mach numbers for thin lip 


0.73 
0.725 
0.72 

Ah 

0.715 
0.71 
0.705 
0.7 

Figure 9: Time history of unsteady nondimensional pres- 
sure at nozzle exit ( x = 2.0, y = 0.642), Mj = 1.15. 


Thick Lip Nozzle 



Figure 10: Sound pressure level at nozzle exit ( x = 
2.0,7/ = 0.642), Mj = 1.15, 326 Hz binwidth. Note 
the strong subharmonics to the left. 
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Figure 1 1 : Frequency of thick lip nozzle over a range of 

Mj. 


and thick lip nozzle geometries. The resulting screech 
tones are analyzed. 

Truly unstructured grids instead of triangulated struc- 
tured grids are used to show the flexibility of the CE/SE 
method and to test its crucial capability of using unstruc- 
tured grids to capture the physics in aeroacoustic prob- 
lems with a feedback loop. The shock cell structure and 
many acoustic aspects of the flows compare well with the 
available experimental data, [4-8] and previous numeri- 
cal results, [15] including the variation of screech fre- 
quency with jet Mach number for the A\ and A 2 modes. 

The thick lip case shows a strong fluid flow entrain- 
ment which produced a counter rotating vortex just up- 
stream of the nozzle lip. This vortex is not seen in the 
thin lip case. 

The difference in the total number of computational 
cells used for a given domain by the solver with a trian- 
gulated structured grid and the truly unstructured grid is 
not significant, so it is possible to consider aeroacoustic 
problems with more complex geometries using truly un- 
structured grids. Great care must be exercised in generat- 
ing an unstructured grid. There were several cases where 
lack of appropriate grid resolution in critical regions of 
the computational domain produced a flow solution but 
no acoustics. 

It is necessary to control the grid clustering to cap- 
ture the aeroacoustics in these jet flows with a feedback 
loop. Grid qualities like smoothness, size, and shape, 
are important for the dissipative aspects of the numerical 
method. 
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Figure 12: Sound pressure level of thick lip nozzle over 
a range of Mj . 
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